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ABSTRACT 



We obtain time-series spectrophotometry observations at the VLT with the aim of partially identifying the dominant oscillation modes 
in the rapidly pulsating subdwarf B star EC 20338-1925 on the basis of monochromatic amplitude and phase variations. 
From the data gathered, we detect four previously known pulsations with periods near 147, 168, 126 and 140 s and amplitudes between 
0.2 and 2.3 % of the star's mean brightness. We also determine the atmospheric parameters of EC 20338-1925 by fitting our non-LTE 
model atmospheres to an averaged combined spectrum. The inferred parameters are T^j = 34,153±94 K, logg=5.966±0.017 and 
logiV(He)/N(H) = -1.642+0.022, where the uncertainty estimates quoted refer to the formal fitting errors. Finally, we calculate the 
observed monochromatic amplitudes and phases for the periodicities extracted using least-squares fitting to the light curves obtained 
for each wavelength bin. 

These observed quantities are then compared to the corresponding theoretical values computed on the basis of dedicated model 
atmosphere codes and also taking into account non-adiabatic effects. We find that the quality of the data is sufficient to identify the 
dominant pulsation at 146.9 s as a radial mode, while two of the lower amplitude periodicities must be low-degree modes with £=0-2. 
This is the first time that monochromatic amplitudes and phases have been used for mode identification in a subdwarf B star, and the 
results are highly encouraging. 



1. Introduction 

Understanding the formation of subdwarf B (sdB) stars is one 
of the remaining challenges in stellar evolution theory today. 
As compact, low-mass, evolved objects located on the Extreme 
Horizontal Branch (EHB) it is generally accepted that they 
descend from the red giant branch, but loose too much of 
their hydrogen-rich envelope to ascend the asymptotic giant 
branch after core He-exhaustion (for a compr ehensive recent 
review on subdwarf B stars see iHeberf [2009). However, the 
details surrounding the mass loss in particular are debatable, 
and a number of evolutionary scenarios have been modelle d 
(e.g. lDorman"eTaT][T99llHan et alJl2Q02l I2003I: I^T&Ij|l2009b . 
These include both single star and binary formation channels 
such as common envelope evolution, stable and unstable Roche 
lobe overflow, and the merger of two Helium white dwarfs. Each 
of these channels should leave a distinct signature on the inter- 
nal composition and the mass distribution of the resulting sdB 
star population. While the latter may be computed dynamically 
for the very rare case of an eclipsing binary, the former can be 
probed only on the basis of asteroseismology. 

Over the last decade, the asteroseismological exploitation of 
subdwarf B stars has become increasingly successful. Full as- 
teroseismological analyses leading to a precise estimate of the 
mass and hydrogen envelope thickness have now been presented 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (proposal ID 
083.D-0415). 



for ten rapidly pulsating sdB stars (see, e.g. Charpinet et alj 
l2008bHvan Grootel et alJl2008tlRandaU et alJl2009i for the most 
recent results). These objects, also known as EC 14 026 or V361 
Hya stars after the prototype ( Kilkenny et aljfl997l) exhibit low- 
order, low-degree p-mode oscillations on a typical time-scale of 
100-200 s, thought to be excited by a classical /(-mechanism 
assoc iated with a loca l overabundance of iron in the driving re- 
gion (ICharpinet et alJ[l99 6. 1997). While the asteroseismologi- 
cal results are beginning to indicate first evolutionary trends (see 
Charpinet et al. 2009, for a review), it is of course highly desir- 
able to verify the accuracy of the claimed structural parameters 
using an independent means. 



One of the most promising ways of doing this is to (partially) 
identify the modes of pulsation observationally. In the currently 
employed "forward method" in asteroseismology, the observed 
periodicities are fit to those computed for a grid of models, and 
an "optimal" model is identified using x 2 minimisation. During 
the period fitting process, the observational frequencies are nat- 
urally associated with a theoretical mode of a given degree index 
I and radial order k (unless there is evidence for relatively fast 
rotation in the target, the computed period spectra are degenerate 
in azimuthal index m). Therefore, any independent knowledge of 
the modal indices can be used to reject an "optimal" model or, 
more usefully, be employed as an a priori constraint for the ex- 
ploration of model parameter space. In certain cases, such con- 
straints may be necessary to isolate just one "opti mal" model 
from a number of very different possible solutions dRandall et al.l 
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2009 ), in others they m ay simply confirm the original solution 
(Ivan Grootel et alJl2008l) . 

The degree index of a mode can be observationally inferred 
using several techniques. For the special case of a rapid ro- 
tator, the presence or absence of rotational splitting can give 
clues to the € value o f at le ast the higher amplitude modes 
(ICharpinet et al.l I2005L 12008a). Other, more universally appli- 
cable methods rely on the wavelength-dependent behaviour of 
the target flux during a pulsation cycle. Several recent projects 
have focussed on the interpretation of line profile variations ob- 
served in EC 14026 stars from high-resolution time-series spec- 
trosco py dTelting et al.l2008t IVuckovic et"aT1l2009t: iTelting et al.1 
2010), but these efforts were always hampered by the low S/N 
of the data. Indeed, the relative faintness of subdwarf B stars 
(V > 12) coupled with the short exposure times acceptable (< 
30 s) for these fast pulsators make this type of study a challenge 
even when granted access to the world's largest telescopes, and 
none of the observations obtained to date were of sufficient qual- 
ity to yield a unique identification of t. Meanwhile, numerous 
studies based on low-resolution time-series spectroscopy of EC 
1402 6 stars (e.g. iJefferv & Pollaccdl2000l lO'Toole et alJ feOOS: 
Telting & 0stensen 2006, to name just a few) have focussed 
mainly on measuring the radial velocity and equivalent width 
variations, and deducing the apparent changes in effective tem- 
perature and surface gravity over a pulsation cycle. Any attempts 
at mode identification were unfortunately inconclusive. 

A more promising approach is to analyse the amplitude and 
phase variations of a pulsation mode as a function of wave- 
length. This has been done quite successfully on the basis of 
multi-colour photometry. While the S/N of the observations is 
still a major limitation, particularly for fainter targets and/or 
low-amplitude pulsations, unique identifica tions of the degree 
index I have been possible in several cases (Jeffery et al. 2004; 
Ran dall et al.ll2005l IChar pinet et al. 2008a; Baran et al. 2008). 
By far the nicest results to date were obtained for the brightest 
EC 14026 star Balloon 090100001 on the basis of UBV light 
curves of extraordinary quality obtained with the three-channel 
photometer LaPoune at the CFHT Not only was the dominant 
pulsation clearly found to be a radial mode, but a further 8 fre- 
quencies could addi tionally be unambiguou sly identified as €=\, 
1=2 or €=A modes (ICharpinet et alJl2008al) . Beyond its intrinsic 
interest, this result proved beyond doubt the validity and accu- 
racy of the theoretical colour-amplitudes calculated according to 
the method outlined in R andall et al.l (120051) . 

The analysis presented here makes use of the same theo- 
retical fram ework and nume rical tools as o ur previous studies 
dRandallet all 120051: ICharpinet et al.ll2008al) . but it is based on 
time-series spectrophotometry rather than multi-colour photom- 
etry. To our knowledge this is the first time that this method of 
mode identification has been applied to a pulsating subdwarf B 
star. While the idea had been present in the back of our minds 
since the development of the the oretical tools, it was the very 
promising results pr esented by Ivan Kerkwrjk et all d2000l) and 
IClemens et al] d2000l) for the ZZ Ceti white dwarf G29-38 that 
finally motivated us to apply for the necessary observing time. 
In that study, just under five hours worth of Keck LRIS time- 
series spectrophotometry were sufficient to measure the degree 
indices for six modes. Of course, G29-38 is an exceptionally 
bright (V ~ 13) white dwarf with several relatively high ampli- 
tude (A > 1%) periodicities, and the monochromatic flux be- 
haviour of g-mode pulsations in ZZ Ceti stars is somewhat dif- 
ferent toJha^_of_£^riodes in EC 14026 stars (cf. Randal l et al.1 
2005; Brassar d et al.ll 1995b . Also, the strong broadening of the 
Balmer lines in white dwarfs greatly facilitates the observational 



Table 1. Time-series spectroscopy obtained for EC 20338-1925 



Night Date (UT) Start time (UT) Length (h) spectra pairs 



1 


2009-08-22 


23:57 


3 


2009-08-25 


00:18 


4 


2009-08-25 


23:48 



04:48 
04:27 
04:52 



654 
615 
656 



characterisation of the flux behaviour across the line profiles, as 
do the longer pulsation periods and correspondingly higher inte- 
gration times acceptable. Thus it was with contained optimism 
that we embarked on the present study. 

Given that the observations envisaged would be obtained 
with the Very Large Telescope (VLT) run by ESO on Cerro 
Paranal, Chile, we selected a target located in the southern hemi- 
sphere. EC 20338-1925 is a moderately bright (V ~ 15.67) 
sdB star with no photometric indication of a cool compan- 
ion ( Kilkenny et al. 2006), and atmospheric parameters around 

7W ~ 35,500 K, log g ~ 5.75, and log/V(He)/N(H) 1.71 

(Ostensen et al.ll2010l) . It was d iscovered to be l ong to the class 
of rapid EC 14026 variables bv lKilkennv et all (120061) . who de- 
tected five periodicities near 168, 151, 147, 141 and 135 s on the 
basis of white-light photometry gathered in 1998. The period 
spectrum was very clearly dominated by the 147-s oscillation, 
which exhibited a high amplitude of nearly 2.5%. This for an 
EC 14026 star unusually strong oscillation mode was the main 
reason we chose EC 20338-1925 over other, brighter target 
candidates. After the observing time had already been granted 
we learned that subsequent observations obtained from 1999 to 
2007 had revealed quite striking amplitude variations for some 
of the modes. In particular, the 147-s oscillation had intermit- 
tently been measured at an amplitude as low as 0.3% (Kilkenny 
2010), while some other frequencies appeared to have constant 
strengths. Such amplitude variations have been observed in a 
number of EC 14026 stars, however their origin is unclear, mak- 
ing EC 20338-1925 an enigmatic target for detailed study. 

2. Observations and Data Reduction 

The observations presented here were obtained using the HIT- 
MS mode of FORS2 mounted on Antu (UT1) at the VLT. This 
rarely used mode is ideal for fast time-series spectroscopy as it 
allows 42 observations per readout, achieved by exposing only 
a small part of the chip and successively shifting the exposed 
blocks of columns across the CCD while already integrating on 
the next spectrum. Moreover, up to two closely spaced targets 
can be observed, allowing for both a main target and a com- 
parison star to be monitored simultaneously. For more detailed 
information on FORS2 and the HIT mode please see the User 
Manual, available online via the ESO webpages. 

We were allocated four consecutive half-nights from 22- 
25 August 2009 for our proposed study of EC 20338-1925. 
Making use of the blue-efficient E2 V CCD formerly installed on 
FORS 1 together with the cross-disperser grism 600B we hoped 
to obtain spectra covering the 3400-6500 A wavelength range. 
Unfortunately, we ran into severe technical problems at the be- 
ginning of the second night of observations, as one of the E2V 
chips of the mosaic died. Since in the HIT-MS mode each spec- 
trum is projected along a block of columns spread across the 
two chips, this resulted in the loss of half the spectral wave- 
length range for the remainder of the observing run. To make 
matters worse, the mask we had cut before the beginning of the 
run placed our target on the dead chip, meaning that we could 
not resume our observations until a new mask had been cut the 
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3250. 3940. 4630. 5320. 6010. 6700. 

Wavelength (A) 

Fig. 1. Typical individual spectra obtained on 22 August 2009 
for EC 20338-1915 (bottom) and the main sequence compari- 
son star (top). The spectra have been arbitrarily normalised for 
visualisation purposes, and are not flux calibrated. 



next day. Thus, despite good atmospheric conditions through- 
out the run we lost the second night completely, and obtained 
only half the wavelength range expected on the third and fourth 
nights. For a log of the observations please see TableQ] 

Since the main aim of our study was to examine the 
monochromatic properties of the target flux as a function of 
time, we avoided slit losses by choosing the largest mask open- 
ings available in HIT-MS mode (5x5"), and essentially obtained 
spectrophotometry as the seeing was always well below the slit 
width. As a natural side effect, the spectral resolution of our data 
is not fixed but seeing-limited, and varies as a function of time. 
The typical seeing on the three nights we obtained data was 
~ 1.1 "at ~4500 A, which with the 600B cross-disperser corre- 
sponds to a wavelength resolution of ~6.5A. On the first night 
(when both chips were working) we covered the 3500-4950 A 
(chip 2) and 5100-6600 A (chip 1) ranges, and opted for the 
bluer 3700-5200 A regime after chip 2 failed. That way, we in- 
cluded the Balmer lines from Ha upwards on night 1, and from 
HjS upwards on nights 3 and 4 for our target star. We also mon- 
itored a F7V comparison star, HD 196286, picked purely for its 
location (7.5' from the target) and brightness (V ~ 10.1, from 
Simbad). The spectral range covered for the comparison star is 
shifted ~400 A bluewards with respect to the target star due to 
its different location on the chip. We set the exposure time for 
each pair of spectra to 25 s, yielding a total integration time of 
1025 s (for 41 measurements) on the CCD before each 36 s read- 
out. This means that only 0.3% of the time-series is taken up by 
dead time in between exposures! 

Given that the HIT-MS mode is not supported by the FORS 
pipeline, we wrote our own IRAF script to batch-process the 
data. This script essentially treats the 82 (41 each for the tar- 
get and comparison star) spectra on each image as echelle or- 
ders, and writes the reduction products out as individual FITS 



files with a time-stamp corresponding to the middle of the ex- 
posure calculated for each pair of spectra. The steps performed 
are bias subtraction, bad pixel interpolation, cosmic ray clean- 
ing, background fitting and subtraction, tracing and extraction of 
the spectra, and finally wavelength calibration. This last step was 
somewhat challenging, since we did not manage to take useful 
arc-lamp calibration data with the available wide-slit masks. In 
the end, the only way to obtain a meaningful wavelength scale 
was to use the spectral lines of the target and comparison star 
themselves, and compute the wavelength solutions for each night 
on the basis of the respective time-averaged combined spectra. 
Luckily, subdwarf B stars show a number of easily identifiable 
Balmer and Helium lines at short wavelengths, and we were thus 
able to obtain a reasonably accurate (within a few A) wavelength 
calibration for at least the blue part of the target spectrum. This is 
sufficient for our purposes, since we do not attempt to detect ra- 
dial velocity variations arising from the pulsations due to the low 
resolution of the spectra and the arbitrary wavelength variations 
induced by spatial shifts of the target across the wide slit. The 
red part of the target spectrum is largely devoid of features, and 
the wavelength solution obtained cannot be trusted for quantita- 
tive analysis. For the comparison star, we were able to identify 
several lines on both chips, and achieve a rough wavelength cal- 
ibration. However, since these measurements are only used in 
wavelength-integrated form in what follows, this calibration is 
needed as nothing more than a guideline. 

Figure Q] shows a typical individual spectrum for both the 
target and comparison star, taken from the first night where both 
chips were available. Note that the comparison star data were 
chopped at the short-wavelength end of the red chip as this 
part of the spectrum was saturated. The high quality is striking 
considering the rela tive faintness of the target (V~15.67 from 
Kilke nny et al.l 2006) and the short exposure time used (25 s). In 
the central part of the blue target spectrum, we measure a signal- 
to-noise S/N ~ 55, while it is lower in the red part at S/N ~ 25. 
For the comparison star an accurate measurement of the S/N is 
complicated by the numerous metal lines, but we can roughly 
estimate S/N ~ 90 for the central part of the blue spectrum and 
S/N > 100 for the red part. It is apparent from looking at the rela- 
tive depth of Ha compared to the other Balmer lines in the target 
spectrum that the red part is a lot less useful than the blue part 
in terms of S/N of the lines, mostly due to the intrinsically low 
flux of sdB stars at longer wavelengths. Adding to this the fact 
that we were able to obtain only half a night's worth of red data, 
and that an accurate wavelength calibration was not possible, 
they become unsuitable for a detailed interpretation. Therefore, 
the analyses presented in the following sections are based on the 
blue data only. 

3. Data Analysis 

3.1. Frequency Determination 

The first step in the data analysis was to obtain pulsation fre- 
quencies for our data set, and at the same time assure the quality 
of the measurements retained for further analysis. In order to 
maximise the S/N we integrated the flux of each spectrum over 
the wavelength range of interest (3650-4950 A) and produced 
broadband fluxes for the target as well as the comparison star. 
These were then used to compute normalised differential light 
curves additionally corrected for differential extinction by fitting 
a second-order polynomial to each night's data. Individual out- 
liers (usually associated with bad seeing) were removed from 
the light curves prior to Fourier analysis, and the correspond- 
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Table 2. Periodicities extracted for EC 20338-1925. We list the 
period and amplitude obtained from our data set, as well as the 
ranges of amplitude measured by D. Kilkenny from 1998-2007. 
The uncertainty on our amplitudes is ~0.05% from the least 
squares fit to the light curve. 



Rank 


Period (s) 


Amp (%) 


A-Kilkemn(%) 


1 


146.940 


2.30 


0.32-2.40 


2 


168.435 


0.45 


0.33-0.46 


3 


126.382 


0.38 


0-0.33 


4 


140.588 


0.20 


0-0.24 



ing spectra were discarded. While we initially included the red 
data, we found that this gave results of lower quality than using 
the blue part of the spectra alone, and thus decided to consider 
only the latter. In addition, the last -hour of data taken on the 
first night had to be disregarded because the stars had shifted 
towards the edge of the slit and some flux loss had occurred, 
making these measurements unsuitable for time-series flux anal- 
ysis. 

We show the Fourier transform (FT) of the combined light 
curve in Fig. [2] By pre-whitening we were able to detect four 
periodicities above the imposed Act threshold for credible pulsa- 
tions (0.2% of the mean brightness), listed in Table [2] together 
with their amplitudes. Three of these, including the strong dom- 
inant pulsation (146.9 s), c orrespond well to those measured in 
1998 (JKilkennv et"ai]|2006l) . while the fourth (126.4 s ) was pre- 
vious ly only detected from a dataset taken in 2003 (Kilkenny 
2010). As briefly mentioned in the Introduction, the amplitudes 
of some of the periodicities observed for EC 20338-1925 (par- 
ticularly the dominant 146. 9-s peak) are highly variable on a 
time-scale of several years, whereas others (such as the 168.4-s 
pulsation) show a remarkably stable amplitude over the course of 
the last decade. This is visualised in Fig. [3] where we have plot- 
ted the amplitudes measured by D. Kilkenny during 7 observing 
runs between 1998 and 2007 for the four oscillations detected 
from our data together with our more recent results (see also the 
last column of Table [2]). It is not clear whether these observed 
amplitude variations are due to real intrinsic amplitude modula- 



m*m)kMMWt 



MwwMLkkkM*', 
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Frequency (mllz) 

Fig. 2. Fourier transform computed on the basis of the integrated 
flux from all blue spectra of sufficient quality gathered during 
the observing run. 



tions, or caused by the beating of very close unresolved modes. 
We hope to shed some light on this in the following sections. 

3.2. Atmospheric analysis 

The second step in the data analysis was to determine the at- 
mospheric parameters of EC 20338-1925, which are necessary 
for the modelling of the monochromatic flux behaviour (see be- 
low). For this purpose, we made a combined spectrum on the 
basis of 995 (blue) spectra spread over the three nights that have 
central wavelength FWHMs close to the average seeing during 
the observing run (between 1.0—1.15")- The result was a time- 
averaged spectrum with wavelength resolution 6.4+0.2 A and 
high S/N of ~400 at the central wavelength. 

We fit this averaged spectrum with two grids of synthetic 
spectra computed from NLTE model atmospheres specially de- 
signed for subd warf B stars using the public codes TLUST Y and 
S YNSPEC (iHubenv & Landll995b lLanz & Hube"nvlll995h . The 
spectra in the first grid contain only H and He in the atmospheric 
composition ("no metals"), while those of the second set ("with 
metals") have a chemical composition made up of H, He, C (0. 1 
solar), N (solar), O (0. 1 solar), S (solar), Si (0. 1 solar) and Fe (so- 
lar). This mixture was derived from the work of Blanchette et al] 
(2008), who determined the abundances of several astrophysi- 
cally important elements for five representative subdwarf B stars 
on the basis of FUSE spectra in conjunction with NLTE model 
atmospheres. The abundance patterns were found to be very sim- 
ilar in all cases, therefore we deemed it appropriate to include the 
most abundant heavy elements in our model atmospheres at the 
level measured. 




1998. aooo. aooa. 2004. aooe. aoos. 2010. 



Fig. 3. Measured amplitudes of the 4 periodicities listed in Table 
|2]over the last 12 years. Apart from the most recent data points, 
which were taken from the observations presented here, the mea- 
surements were kindly made available to us by D. Kilkenny. The 
errors were not provided, but should be similar to or smaller than 
those shown for our data. Zero amplitudes imply that the pulsa- 
tion was not detected in a particular data set down to a threshold 
of -0.1%. 
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Fig. 5. Left panel: Available H and He lines from the combined FORS2 spectrum of EC 20338-1925 (thin lines), overplotted with 
the best atmosphere fit not including metals (thick lines). The atmospheric parameters derived are indicated in the plot. There are 
signs of interstellar absorption as revealed by the Ca II K feature in the blue wing of He. Right panel: The same as the left panel, 
except that here metals are included in the proportions listed in the text. 
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Fig. 4. Normalised averaged spectrum for EC 20338-1925 (thin 
line) overplotted with the best synthetic spectrum "with metals" 
(thick line). The vertical dashed lines indicate the positions of 
the SHI, Oil and Nil lines thought to be present in the observed 
spectrum. The atmospheric parameters derived from this fit are 
also given. 



Fig.|5]shows the fits to the Balmer and Helium lines obtained 
for both the "no metals" and "with metals" cases, and also gives 
the inferred values for the effective temperature, logarithmic sur- 
face gravity and logarithmic fractional Helium abundance. The 
uncertainty estimates given of course refer only to the formal fit- 
ting errors and most certainly underestimate the real errors. As 
is typically the case for EC 14026 stars with T eS -35,000 K, the 
resulting atmospheric parameters are only marginally influenced 
by the inclusion of metal s in the model at mospher es, and ar e 
consistent with each other dHeber et al. 2000; Geier e t al.ll2007h . 
And indeed, the line profile fits shown in Fig. [5] hardly change 
from one plot to the next. 

There are nevertheless some indications that the solution in- 
corporating metals is to be preferred. Despite its relatively low 



resolution, the averaged spectrum of EC 20338-1925 does hint 
at the presence of a few heavy elements, thanks to the high qual- 
ity of the data. In particular, the normalised spectrum plotted 
in its entirety in Fig. [4] contains lines corresponding to absorp- 
tion from Nil £4237.05 and 4241.79 A), Oil (4253.87 A) and 
SHI (4284.97 A). These are marked by vertical dotted lines in 
the plot, and are recovered quite nicely by the overplotted "with 
metals" synthetic spectrum. Another indicator for the presence 
of metals in EC 20338-1925 is the relatively poor fit to the He 
II line at 4686 A (see Fig. [5}. While the model incorporating 
metals fares better in the matching of this line than the syn- 
thetic spectrum without metals, it appears that even the "with 
metals" chemical composition does not include enough metals 
to accurately represent EC 20338-1925. Given the fact that the 
spectroscopic solution is not very sensitive to the metallicity of 
the model atmosphere for our target, this does not constitute a 
problem for the accuracy of the inferred atmospheric parameters. 
Therefore, we adopt the "with metals" values of r e ff=34,153 K 
and log g=5.966 in what follows. 



3.3. Observed monochromatic amplitudes and phases 

The final, and arguably most interesting part of the data anal- 
ysis was the computation of the observed monochromatic am- 
plitudes and phases for the periodicities extracted. For this ex- 
ercise we used only the blue target spectra also employed for 
the computation of the broadband Fourier transform shown in 
Fig. [2] thereby ensuring a high data quality and the exact same 
sampling. The latter is important since we used the periods de- 
termined in Section 3.1 as input for the least-squares fitting pro- 
cedure detailed below. All selected target spectra were chopped 
to the 3650-4950 A range of interest and binned into common 
3.25 A wavelength units, which were chosen to adequately sam- 
ple the average ~6.5 A resolution of the original data. 

Light curves were then obtained for each of the 401 wave- 
length bins, and corrected for differential extinction using the 
comparison star data. In order to achieve the best possible S/N, 
we integrated the comparison star spectra from 3650 A to their 



S.K. Randall et al.: Monochromatic amplitude and phase variations for EC 20338-1925 




3650. 3850. 4050. 4250. 4450. 
Wavelength (A) 



4650. 4850. 




Wavelength (A) 



Fig. 7. Left panel: We show the observed amplitude as a function of wavelength for the dominant pulsation mode /j . For reference 
we plot the averaged spectrum of EC 20338-1925 on an arbitrary y-axis scale and also indicate the central wavelengths of the 
Balmer and prominent He lines (dotted vertical lines). The thicker yellow bars mark the wavelength bins for which the Fourier 
Transforms plotted in Fig.|6]were computed. Right panel: The same as the left panel, but for the observed phase. 
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upper limit of ~ 4700 A and employed the resulting broadband 
fluxes for the calculation of the differential light curves. For each 
wavelength bin, and using the periods listed in Table [2] as fixed 
input parameters, the corresponding amplitudes and phases were 
determined from a least-squares fit to the light curve assuming 
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a sinusoidal signal. While this did not involve the computation 
of Fourier spectra, we nevertheless calculated them for a small 
number of selected wavelength points as a sanity check. For il- 
lustration purposes we show some of these Fourier transforms in 
Fig. [6] Note that here, the spectra were binned into units of 10 A 
to somewhat decrease the noise level. 

Fig.|7]shows the monochromatic amplitudes and phases de- 
rived for the dominant mode f\. For reference we also show 
the averaged spectrum of EC 20338-1925, and indicate the rest 
wavelengths of the Balmer and most prominent Helium lines. 
Looking at the amplitude spectrum we find that the amplitude of 
pulsation in the continuum drops quite rapidly between 3650 and 
4050 A and then stabilises somewhat. This was to be expected, 
since t he same effect is observed from multi-colour photome- 
try (e g. lCharpinet et alj2008bHRandall et al.l2005l:ljefferv et al.1 
2004). In the Balmer line cores, the pulsational amplitude is sig- 
nificantly higher than in the continuum, and again decreases with 
increasing wavelength of the line. Considering the Helium lines 
strong enough to leave a visible imprint on the amplitude spec- 
trum, it appears that the neutral helium lines at 4026 and 447 1 A 
are associated with higher amplitudes compared to the contin- 
uum, while the pulsational signature is extremely weak in the 
He II line at 4686 A. The observed phase on the other hand ap- 
pears to be constant in the continuum, jumps occuring only in the 
line centers. This is consistent with the very small phase shifts 
observed on the basis of multi-colour photometry dJefferv et al.l 
l200Hl2005l) . 

4. Modelling the observed monochromatic 
amplitude and phase variations 

4. 1. Computation of theoretical amplitudes and phases 

For the computation of the theoretical monochromatic ampli- 
tudes and phases we closely follow the approach developed by 
Randall et al. (2005) for the application to multi-colour photom- 
etry. We refer the interested reader to that paper for details on 
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the underlying theory and the models, and only give a very brief 
overview of the methodology followed here. 

Our theoretical amplitudes and phases rely primarily on 
monochromatic quantities derived from model atmospheres. The 
computation of these quantities involves not only the standard 
specific intensities, but also the corresponding derivatives with 
respect to effective temperature and surface gravity across the 
visible disk of the star. They are calculated on the basis of a grid 
of LTE model atmospheres characterised by a uniform composi- 
tion specified by log |W(He)/N(H)] = -2.0, i.e. without metals. 
Whereas it would be interesting to assess the effect that incor- 
porating metals would have on this type of calculation in the 
future, the current grid is quite sufficient for the purposes of this 
study. In any case, the heavy element abundance pattern of EC 
20338-1925 ha s not yet been characterised, and most likely de- 
viates from the iBlanchette et all (2008) composition employed 
for the atmospheric analysis detailed in Section 3.2. We there- 
fore compute the monochromatic atmospheric quantities using 
the existing H/He model grid, assuming an effective temperature 
J e ff=34,153 K and a logarithmic surface gravity logg=5.966 as 
derived above. 

While the monochromatic atmospheric quantities thus com- 
puted are the main contributors to the relative first order per- 
turbations of the emergent flux (from which the monochromatic 
amplitudes and phases directly follow), they neglect to take into 
account non-adiabatic effects. The quantities describing the de- 
parture from adiabacity in amplitude and phase must instead 
be derived from full non-adiabatic pulsation calculations on the 
basis of stellar models. For this purpose, we emplo y our well- 
known "sec ond-generation" models (see, e.g. Charpinet et al. 
1996, 2001) in conjunction with adiabatic and non-adiabatic pul- 



sation codes based on finite-elem ent techniques (Brassard et alJ 
1 19921: [Fontaine & Brassard! fl9 94). These are exactly the same 
tools that have been used very successfully in the past for the 
asteroseismological analysis of EC 14026 stars. They require 
four structural input parameters: T e s , log g, the total stellar mass 
M„ and the fractional thickness of the hydrogen-rich envelope 
logg(H) = logM(H)/M*. In our computations we set the first 
two quantities to those inferred in Section 3.2, while for the 
other two we selected representative values of M* = 0.48 M Q 
and logg(H) = -4.3. The two latter parameters were chosen to 
be typical of rapid sdB pulsators as found from asteroseismolog- 
ical studies, but do not significantly affect the results as long as 
they take on physically reasonable values. 

The quantities describing the departure from adiabacity were 
first computed for the period spectrum of the representative EC 
20338-1925 model, and subsequently determined for the dom- 
inant 146. 9-s dominant mode using interpolation techniques. 
While for p-mode pulsators these quantities are quite sensi- 
tive to the period of the oscillation in question, they are luck- 
ily not dependent on the degree index of the mode, and can 
thus be inferred quite easily for any period in the excited fre- 
quency range predicted by the model. For reference, at the pe- 
riod of the main mode the non-adiabatic coefficients are R = 
0.810 and TV = 2.761 rad compared to their adiabatic values of 
R — 1 and ¥7- — n (for a definition of the parameters please see 
Randall et al. 2005). With both the monochromatic atmospheric 
quantities and the non-adiabaticity coefficients available, we can 
now easily compute the monochromatic amplitudes and phases 
for EC 20338-1925. 

Fig. [8] shows the relative amplitudes as expected for the 
146.9-s mode of EC 20338-1925 for degree indices from ^=0 
to {=5 in the wavelength range of interest. The curves have been 
normalised to unity at 3650 A and were degraded in wavelength 
resolution to the 6.5 A typical of the observational data. Higher 
order modes are not illustrated as they are thought to have ex- 
tremely low amplitudes in the continuum due to cancellation 
effects when integrating across the visible disk o f the star, and 
should thus not be observable. As explained in Randall et al. 
(2005), the theoretical amplitudes and phases are, to first order, 
independent of the azimuthal index m, and can therefore be used 
for resolved or unresolved frequency components of multiplets 
split by slow rotation. The curves quite clearly show that the 
discrimination between low-order modes with I < 2, and par- 
ticularly between radial and dipole modes is challenging due to 
their very similar behaviour as a function of wavelength. This 
has always been the main difficulty also with the interpretation 
of multi-colour photometry and line-profile variability, and has 
limited many studies to the ambiguous identification of the main 
pulsations as low-order (i.e. ( — 0, 1 or 2) modes due to the in- 
suffic i ent S/N of the data (e.g .lJeffery et alj2005H T remblav et al.l 
2006; Vuckov ic et al.l 120 09). Given that the amplitude depen- 
dence is greatest at short wavelengths, it is therefore imperative 
to observe as far in the blue as possible. 

Comparing the curves displayed here for EC 20338-1925 
to equivalent data for other EC 14026 stars (see in particular 
Fig. 2 of iRandalle t al. 2005), it turns out that we were quite 
lucky with our target in terms of discriminative power between 
the amplitu des of low-order mod es. Indeed, as can be seen from 
Fig. 28a of Randall et al.l (120051) . the amplitude-wavelength be- 
haviour of radial, dipole, and quadrupole modes diverges more 
strongly for stars with a higher surface gravity. By contrast, the 
(-3 mode shows a shallower amplitude gradient as a function of 
wavelength. EC 20338-1925 being one of the EC 14026 stars 
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Fig. 10. Left panel: We show the normalised observed amplitude for f\ (black curve) together with the errors on the data points (thin 
curves shifted downwards along the y axis for clarity). The theoretical monochromatic amplitude as calculated for a radial mode 
with P=146.9 s and model parameters characteristic of EC 20338-1925 is overplotted (red curve). Right panel: Similar to the left 
panel, but referring to the monochromatic phase. The theoretical monochromatic phase was computed using the same parameters 
as for the amplitude, and additionally a radial velocity shift of 1 1 km/s was assumed. 



with the highest surface gravity measured is then definitely to 
our advantage, and increases the chances for the unique mode 
identification attempted below. 

4.2. Fitting the observed monochromatic amplitudes and 
phases 

We now compare the theoretical monochromatic amplitudes 
computed in the previous section to the observations. Fig. [9] 
shows the curves displayed in Fig.|8]overplotted on the observed 
monochromatic amplitude of the dominant f\ mode, which has 
also been normalised to unity at 3650 A. From visual inspection 
alone it is obvious that the I — solution presents by far the best 
match to the observational data, followed by the curves for I - 
1 and, with a significant degradation t — 2. The amplitudes for 
the higher degree modes match the data very poorly, and can be 
excluded from the outset. Computing the quality-of-fit Q on the 
basis of x 1 fitting confirms this result, with (2=0.755 for ( = 0, 
Q = 0.340 for € = 1 and Q = 0.056 for 1 = 2. For the higher de- 
gree modes, Q <k 0.001, which i mplies that they c an formally be 
excluded as possible solutions dPress et al.lll986l) . The Q value 
determined for the radial mode implies a very good fit, more 
than twice as good as that for the dipole mode. While the un- 
certainties on the observed amplitudes are too large to formally 
exclude the solutions for I — 1 and I — 2 as possible matches 
to the data, we believe that the far superior quality of the fit for 
I — leaves no doubt that the dominant pulsation observed in 
EC 20338-1925 is a radial mode. 

In order to illustrate the quality of the fit more clearly, we 
zoom in on the monochromatic amplitude plot for I — in the 
left panel of Fig. [10] Here, we also include the uncertainty on 
the observed amplitudes as obtained from the least-squares fits 
to the light curve for each wavelength bin (thin vertical line seg- 
ments shifted downwards along the y-axis for clarity). It can be 
seen that, while the observed amplitude behaviour in the contin- 



uum is matched nearly perfectly by the theoretical curve, there 
are some discrepancies in the spectral lines. Most strikingly, the 
amplitude dip corresponding to the He II line at 4686 A is pre- 
dicted far weaker than observed, much like what we found for 
the model atmosphere fit to the time-averaged spectrum in Fig. 
[5] This is almost certainly due to the lack of metals in the model 
atmospheres used for the computation. A more puzzling feature 
is that the theoretical amplitudes in the Balmer line cores system- 
atically deviate from those observed: while they are similar to, 
or slightly lower than the measurements for the higher Balmer 
lines (above He), they increasingly overestimate the observed 
amplitudes in the lines at higher wavelengths (for He, H5, and 
especially Hy and H/3). We are not sure whether this is due to 
some observational effect or caused by missing ingredients in 
our models (particularly concerning the chemical composition 
of the atmosphere). At the request of the referee we re-computed 
the g-value of the fit taking into account only the continuum 
and the extended line wings (specifically, we masked out bands 
±6Awide centered around HyS, Hell 4686, Hy, US, He and H8). 
As expected, the quality-of-fit is improved, with 2=0.999 (^=0), 
2=0.503 ((=1), 2=0.223 (f=2) and Q «0.001 for higher lval- 
ues. This is simply associated with the fact that there are now 
less deviations between the predicted and observed amplitudes, 
and nicely illustrates the very good match achieved for a radial 
mode in the continuum. 

While we are conviced that the dominant mode in EC 
20338-1925 is radial from the amplitude data alone, we also 
analysed the corresponding phase data as a consistency check. 
Given the high noise level of the measurements, we assumed a 
degree of I — from the outset, and attempted to fit the theoret- 
ical phases to those observed. In order to be able to compute 
the monochromatic phases accurately, we extended our exist- 
ing code for the computation of the pulsational flux perturba- 
tion by an optional radial velocity (RV) parameter. This parame- 
ter was not necessary for the theoretical quantities presented by 
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iRandall et al.l (|2005) because in that study we were interested 
only in wavelength-integrated phase measurements, which are 
dominated by the contributions from the continuum. As can be 
seen from Fig. [7] these yield a relatively constant phase as a 
function of wavelength, which is consistent with the very small 
phase shifts predicte d from one filter bandpass to the next in 
IRandall et al.l (120051) . However, in a spectral line the moving 
matter produces significant distortions and phase delays over 
a pulsation cycle, and these need to be taken into account for 
monochromatic phase calculations. The distortions also affect 
the amplitudes in the spectral lines, but the effect is thought to be 
relatively small for the broader lines and was not included in the 
computations presented above for technical reasons. It remains 
to be investigated whether an incorporation of the Doppler shifts 
in the monochromatic amplitude computations could help im- 
prove the match between observations and theory in the Balmer 
lines. 

The RV parameter incorporated into the phase calculations 
corresponds to the measurement of the semi-amplitude of the 
RV integrated over the disk, the same quantity that has been 
the subject of numerous time-series spectroscopy studies to 
date. Unfortunately, our data are simply not suitable for such 
an analysis due to the relatively low resolution, problems with 
the wavelength calibration (see Section 2), and spurious wave- 
length shifts caused by the star moving within the rather wide 
slit. We therefore have no observational constraint on the RV 
parameter, and instead treat it as a free parameter while fit- 
ting the observed monochromatic phases. That way, we were 
able to model the observed phases quite nicely, as can be 
seen in the right-hand panel of Fig. [10] Note that the y-axis 
now measures the phase in radians rather than in seconds as 
in Fig. [7] for easier comparison to the computed values. This 



causes an apparent inversion of the curve due to a minus sign 
in the relation between phases specified in radians and those 
quoted in seconds. The radial velocity parameter used in the 
plot has a semi-amplitude of 11 km/s, which is quite simi- 
lar to the values found for other high-amplitude EC 14026 
pulsations from RV studies (e.g. ~ 12 km/s for the dominant 
mode in PG 1325 + 101 from lTelting & 0stensenll2OO4 19 km/s 
for the un usually strong pulsation in Balloon 09010001 from 
Telting & 0stensen 200 6J and 15 km/s fo r the main frequency 
in PG 1605+072 from lO'Toole et alJl2005l) . We believe that this 
consolidates our interpretation of the dominant pulsation as a ra- 
dial mode. Nevertheless, it is clear that given the observational 
limitations on the S/N as well as the lack of an independent ra- 
dial velocity estimate, the discriminative power for the degree 
index of the mode lies in the amplitude, rather than the phase 
information. 

Considering the lower amplitude modes extracted from our 
data, the monochromatic quantities are a lot more noisy, mak- 
ing a unique mode identification impossible. Nevertheless, we 
repeated the exercise of overplotting the normalised theoretical 
amplitude curves on the measurements for fa and fo (the data for 
ft, were too noisy due to the much lower amplitude). We used 
exactly the same theoretical values as for the dominant mode, 
because the small variation due to the difference in period is 
negligible when dealing with data this noisy. The results of the 
overplotting are shown in Fig.QT] Unsurprisingly, the S/N of the 
observations is too low to distinguish between the low-order de- 
gree modes with ( — 0, 1 and 2, however the higher degree modes 
present an inferior fit, and can quite safely be excluded. The cor- 
responding monochromatic phases meanwhile are too noisy to 
gain any information whatsoever. 



10 



S.K. Randall et al.: Monochromatic amplitude and phase variations for EC 20338-1925 



1.5 



1.0 

o 

" 0.5 



> 1.5 



U 1.0 




0J 

2 



a 



<0 1.5 

oj 

g 1.0 

O 

0.5 







Wavelength (A) 

Fig. 9. Theoretical monochromatic amplitudes for modes with 
I - 0-5 overplotted on the observed amplitudes for the domi- 
nant f\ mode in EC 20338-1925. The curves have all been nor- 
malised to unity at 3650 A. 

5. Conclusion 

In this study we attempted a mode identification of the main pul- 
sation frequencies in the rapid sdB pulsator EC 20338-1925 on 
the basis of monochromatic amplitude and phase variations. We 
were able to obtain 3 half-nights of time-series spectrophotom- 
etry with the VET instrument FORS2 in HIT-MS mode and de- 
tect four pulsations above the 4-cr threshold. The S/N of the data 
was sufficient to extract high-quality monochromatic amplitudes 
and phases only for the dominant pulsation, noisier results be- 
ing obtained for the remaining periodicities. Using full model 
atmosphere codes appropriate for the atmospheric parameters of 
our target and also incorporating non-adiabatic effects, we com- 
puted accurate theoretical amplitudes and phases for comparison 
to the observations. From the fits achieved to the observed am- 
plitudes alone we determined the dominant pulsation to be a ra- 
dial mode. This conclusion was found to be consistent with the 
monochromatic phase shifts observed under the assumption of 
a very reasonable radial velocity variation of 1 1 km/s. For two 
lower amplitude periodicities the data were found to be of suffi- 
cient quality to exclude higher degree modes with t > 3. 

To our knowledge, this is the first time that mode identifica- 
tion has been attempted for an EC 14026 star using monochro- 
matic amplitudes. The fact that we were able to identify the 
main periodicity as a radial mode is therefore quite encourag- 
ing for future studies, particularly considering the relative faint- 
ness of our target. Moreover, the excellent match between theo- 
retical predictions and the observations in the continuum con- 
firms the basic validity of our approach and our models, al- 



though the discrepancies in the line cores indicate that there is 
still room for improvement, especially with respect to the chem- 
ical composition assumed in the model atmosphere calculations. 
Unfortunately, sdB stars are all chemically peculiar and con- 
tain heavy elements in varying abundances, implying that for 
a perfect match each target would first have to be submitted to 
a detailed abundance analysis and then analysed using a spe- 
cially computed grid of model atmospheres. This is extremely 
time consuming both on the observational and the computational 
front, and not necessary for the practical application presented 
here, given that we have proved mode identification to be possi- 
ble on the basis of the continuum behaviour alone. 

It is quite striking that in the few cases where observa- 
tional mode identification has been possible for rapidly pulsat- 
ing subdwarf B stars, the dominant mode of pulsation was al- 
ways found to be a rad ial mode (for Balloon 090100001 from 



Charpinet et al. 2008b, KPD 2109+4401 from Randall et al. 



2005, HS 2201+2610 from Silvotti et al. 2010, submitted, and 



now for EC 20338-1925). Of course, this cannot automati- 
cally be assumed to hold true for all the highest amplitude fre- 
quencies detected in EC 14026 stars, but it does seem reason- 
able to deduce that mode visibility across the disk of the star 
is a very important factor influencing the observed amplitude. 
Asteroseismological solutions assigning high-degree modes to 
the strongest frequencies observed can in all likelyhood be dis- 
regarded from the outset, as has occasionally been done in the 
past t o distinguish betwe en several possible "optimal models" 
(e.g. Randall et al. 2009). This approach allows unambiguous 
model solutions to be found on the basis of fewer observed pe- 
riods than if all observed frequencies are allowed to take on de- 
gree indices of ^=0,1,2, and 4, and can therefore save valuable 
observing time. 

One puzzling implication of our identification of the 146. 9-s 
periodicity as a radial mode is that this frequency cannot then 
contain closely spaced components due to a rotationally split 
multiplet. Consequently, it is difficult to explain the strong am- 
plitude variations observed over several years (see Fig. 3) in 
terms of the beating of unresolved frequencies. Given that the 
photometry data obtained by Dave Kilkenny have time base- 
lines of several days, the two independent harmonic frequen- 
cies would have to be spaced less than ~ 2 juHz apart in order 
to not be resolved. Additionally, the strong amplitude variations 
observed can be produced only if the two pulsations have com- 
parable amplitudes, effectively implying two very closely spaced 
modes with degree indices t — and/or 1 . From our stellar mod- 
els we know that the required proximity in frequency is impossi- 
ble for modes of the same degree index, and ext remely unl ikely 
for the combination of t — and I = 1 (see, e.g. lCharpinet et atl 
2000). This is particularly true for a target with as high a surface 
gravity as EC 20338-1925, since the period density predicted 
is lower than for a more typical, less compact EC 14026 star. 
Therefore, we believe there is a high probability that the ob- 
served amplitude variations for the main mode are intrinsic to 
the star. 

In conclusion, we have shown that the monochromatic am- 
plitude and phase variations can be effectively used for mode 
identification in rapidly pulsating subdwarf B stars. The most 
critical factor is the quality of the observational data, as the 
S/N level attained limits the discriminative power especially be- 
tween modes of low degree indices with ( < 2. Even though we 
were allocated observing time using a highly suitable instrument 
mounted on one of the world's largest telescopes we were able to 
unambiguously identify the degree index for only the dominant, 
rather high amplitude mode. On the other hand, EC 20338-1925 
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is relatively faint, and there are several brighter targets available. 
Moreover, our observing run was affected by severe technical 
problems, due to which we lost a quarter of our allocated time 
completely, and did not obtain enough red data to constructively 
include it in our analysis. Therefore it seems likely that better 
results could be obtained in future similar studies. 

It is not straightforward to objectively assess the relative 
merits of the different techniques currently in use for mode iden- 
tification in subdwarf B stars, and opinions on this will vary. 
We believe that for the time being the exploitation of the pulsa- 
tional amplitude's signature on I as a function of wavelength 
is the most promising method, and it is so far the only one 
to have yielded unambiguous determinations of the degree in- 
dex t. Given that the work presented here constitutes the first 
study employing monochromatic rather than broadband ampli- 
tudes, it is too early to say which of these two closely related 
techniques is more efficient. One drawback of the amplitude- 
wavelength approach to mode identification is that medium- 
sized or large telescopes coupled with specialised instruments 
are needed in order to obtain rapid time-series data of sufficient 
quality. However, this is equally true for the study of line-profile 
variations, where the low S/N achievable even with the world's 
largest telescopes has so far prevented unambiguous mode iden- 
tification. Radial velocity surveys based on data from smaller 
telescopes have been successfully used to confirm periodici- 
ties previously detected from photometry, but we lack the the- 
oretical framework necessary to use the information gained for 
mode identification. When supplementing simultaneous multi- 
colour photometry, radial velocity measurements can indeed im- 
prove the di scriminative power of the former, as was shown by 
iBaran et a l. (2008, 2010). However, we would like to point out 
that for Balloon 09010001, the only case where a direct com- 
parison of the two techniques is possible, four nights of stan- 
dalone high S/N m ulti-colour pho tometry from a medium-size 
telescope (Charpin et et al.l|200 8b) yielded an identification for 
significantly more modes than the combined results from 30 
hours of simultaneous time-series spectroscopy obtained on a 
medium-size telescope an d 120 hours of mu lti-colour photome- 
try from a small telescope (IBaran et al.ll2008l) . Therefore, we are 
confident that the exploitation of multi- wavelength amplitude in- 
formation along the lines presented here makes the most efficient 
use of telescope time (at least until it becomes possible to obtain 
accurate monochromatic amplitudes and radial velocities with 
a single instrument), and will prove highly conducive to mode 
identification and, consequently, asteroseismology in the future. 
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